arXiv:1509.03635v2 [astro-ph.EP] 16 Sep 2015 


Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 17 September 2015 (MN style file v2.2) 


Triggered fragmentation in self-gravitating discs: forming 
fragments at small radii 

Farzana 

^Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge, CBS OH A, United Kingdom 
Institut fiir Astronomie, ETH Ziirieh, Wolfgang-Pauli-Strasse 21, 8093 Ziirieh, Switzerland 
^ Institut fiir Astronomie und Astrophysik, Universitdt Tubingen, Auf der Morgenstelle 10, 72076 Tubingen, Germany 
^Sehool of Physies, University of Exeter, Stoeker Road, Exeter, EX4 4QL 


17 September 2015 


ABSTRACT 

We carry out three dimensional radiation hydrodynamical simulations of gravitation¬ 
ally unstable discs to explore the movement of mass in a disc following its initial 
fragmentation. We find that the radial velocity of the gas in some parts of the disc 
increases by up to a factor of ~ 10 after the disc fragments, compared to before. While 
the movement of mass occurs in both the inward and outward directions, the inwards 
movement can cause the inner spirals of a self-gravitating disc to become sufficiently 
dense such that they can potentially fragment. This suggests that the dynamical be¬ 
haviour of fragmented discs may cause subsequent fragmentation to occur at smaller 
radii than initially expected, but only after an initial fragment has formed in the outer 
disc. 

Key words: accretion, accretion discs - protoplanetary discs - planets and satellites: 
formation - gravitation - instabilities - hydrodynamics 


1 INTRODUCTION 


Self-gravitating discs are known to be important on a wide 
range of astrophysical scales. If they fragment, a number of 
objects may potenti ally form out of discs o f different sizes 
inclu ding stars (e.g. Navakshin et al.ll2007l b br own dwarfs 
fe.g. IStamatellos et al.ll2007l l and planets (e.g. IBos 
iMaver et, al.ll2nn!ib 

The stability of a self- gravitating disc can be describ ed 
by the stability parameter (| Safronov! 1 1960l : lToomrelll964l b 




ttSC’ 


( 1 ) 


where Cs is the sound speed in the disc, /^ep is the epicyclic 
frequency, which for Keplerian discs is approximately equal 
to the angular frequency, Q, E is the surface mass density 
and G is the gravitational constant. For an infinitesimally 
thin disc to fragment, the sta bility parameter must be less 
tha n a critical value, Q crit ~ 1 (|Safronovlll960l : iToomrell 19641 : 
see iHelled et al.ll20l3 for a recent review). 

Investigations into gravitational instability typically 
consider discs up until the fragmentation stage. Those sim¬ 
ulations that do consider the subsequent evolution of such 
discs focus on the properties of the fragmen ts that form 
and th e longer te r m sta t e of the system (e.g. M^eret^^y 


I 2 OO 2 I : iRice et ^ 120031 : IStamatellos fc WhitworthI 1200' 
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IStamatellos et aP 1201 ll : IZhu et ahl 12012l : IVorobyo^ l2013l f . 
Following fragmentation there are a number of processes 
that can take place. These may affect the fragment (e.g. col¬ 
lapse, accretion, me tal enrichment, tidal disruption, migra¬ 


tion k, gap-opening: iGalvagni et ak 20121: iBolevet al.ll2011 
20R; iNavakshinl 1201 Ol : iBaruteau et al.l 1201 ll: iMalik et a! 


20^ ) , the star (e.g. FU Orion i s outb ursts; Ibodato Clarke 


.2004 ; I Vorobyov BasiJ l2005l . l2006h and the disc (e.g. loss 
of sp i ral structure, gap-opening: | Stamatellos V WhitworthI 
l2009l : lMerull2010l : iMalik et al.ll2015l b However, little is un¬ 
derstood specifically about the movement of the mass in 
the disc following the formation of the first fragment and 
how this affects the possibility of su bsequent fragmenta¬ 
tion. I n the context of planet formation lArmitage V HansenI 
(|l999l f performed Smoothed Particle Hydrodynamics simu¬ 
lations of self-gravitating discs and inserted a planet (using 
a smoothed point mass) to explore the subsequent effect on 
the disc. They found that the presence of a planet triggers 
the formation of fragments both exterior and interior to it. 
These simulations, however, considered an isothermal equa¬ 
tion of state, which favours fragmentation. 


In this paper we use 3D SPH simulations with radia¬ 
tive transfer to self-consistently model the fragmentation of 
a young self-gravitating disc and the subsequent movement 
of gas in the disc. We show that the formation of the first 
fragment causes the radial velocity of the disc mass to in¬ 
crease (by up to a factor of 10) in both the inwards and 
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outwards directions. The inwards movement, which occurs 
because of the presence of the first fragment, can cause the 
inner regions of a disc that are either stable or marginally 
stable (so that Q ~ 1) to become sufficiently dense and un¬ 
stable such that they may potentially fragment. 

In Section we describe our numerical method and 
present the details of the simulations in Section [S] We 
present our results in Section [H and discuss and make con¬ 
clusions in Sections [5] and [6l respectively. 


2 NUMERICAL METHOD 



-100 -50 0 50 100 
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The Smoothed Particle Hydrodynamics (sph) code 
used to carry out the simulations presented here 
is almost exa ctly the same as the code use d by 


originally developed bv iBenz 


iMeru Bate! (|2010l ) 

S and further developed bv [ Bate. Bonnell Sz Price 
_2 _ Whitehouse. Bate Mon aghanI (|2005h . 


IWhitehouse Bat3 (j2006l l and IPrice Batel (|200'^ b We 
use the flux-limited diffusion approximation to model the 
radiation transport with two temperatures (that of the gas 
and the radiation field). To m odel the shocks in the discs , 
we use an artificial viscosity ([Monaghan k, Gingoldl Il983l i 
with SPH parameters fixed at (<asPH,/^ sph) = (0.1,0.2). 
W e use grey Rosseland mean opacit ie s bas ed on the values 
of IPollack. McKay Christ offerson| ( 198^) for interstellar 
molecular dust grains and lAlexand^ ( 1975li for the gas con¬ 
tribu tions at high temperatures (see IWhitehouse Batel 
I 2 OO 6 I for details on the opacity tables used). We carry 
out simulations with opacity values that are lower than 
interstellar values and scale t he values in t h e tab les by the 
required factor as done so bv iMeru Batel (|2010l b 

The equation of state used in these simulations assumes 
that the gas is composed of hydrogen (70%) and helium 
(28%) and includes the rotational and vibrational modes of 
molecular hydrogen, dissociation of molecular hydrogen and 
ionisation of both hydrogen and helium. We use a ratio of 
specific heats, 7 = 5/3. Our models assume that the pres¬ 
ence of dust only affects the value of the opacity in the disc. 

We model the effects o f stellar irradiation u sing a two- 
layer approach, as done bv iMeru Batel (|201Ql b The mid¬ 
plane region is modelled using the flux-limited diffusion ap¬ 
proximation and the boundary layer of particles maintains a 
fixed temperature profile (assumed to be set by the radiation 
field of the star) such that any energy that is passed to these 
boundary particles is effectively radiated away. The vertical 
location of the boundary between the optically thick part 
of the disc and the o ptically thin atmosph ere is determined 
slightly differently to iMeru Batel (|201Ql i : their boundary 
location was determined at the start of the simulation and 
then kept fixed throughout while our method re-evaluates 
the boundary location at least 33 times per outer rotation 
period (ORP). The boundary height is the maximum out 
of: the height above the midplane where the optical depth, 
T = 1, or the height above the midplane where 10% of the 
particles at any radius are in the boundary. Therefore, the 
boundary is at the vertical position where the optical depth, 
T < 1, and our choice ensures that the “bulk” of the disc (i.e. 
consisting of particles that lie closer to the disc midplane) 
is simulated using radiative transfer, r ather than simplified 
energetic calculations. As described bv iMeru Batel |2010l ) 


Figure 1. Logarithm of the gas temperature (in K) rendered 
in cross-sectional views of the simulations before any fragments 
form. Top panel: Simulation 1 at t = 1.29 ORPs. Bottom panel: 
Simulation 2 at t = 1.92 ORPs. Axis units are in an. The surface 
temperatures of the discs are 10 — 50 K, and the outer disc is 
close to being vertically isothermal. 


the particles that make up the boundary do change over 
time since particles may move across the boundary from the 
optically thick to optically thin region, or vice versa. 

To decrease the computa tional expense, we turn any 
fragments into sink particles ([Bate. Bonnell Sz Pried [T995h 
when they reach a critical density of pcrit = 7 — 9 x 
10“®g/cm^, i.e. denser than that required for second core 
collapse, allowing us to follow the evolution of the disc fur¬ 
ther. In the simulations presented here, the first fragments 
form at ~ 33 — 55 AU and have masses ~ IMj. This gives a 
Hill radius (defining the region from which a core typically 
accretes from) of ~ 2 — 3AU. We therefore model the frag¬ 
ments using an accretion radius of 1 AU i.e. we assume that 
any gas within 1 AU of the sink particle can be accreted 
onto it. 


3 SIMULATIONS 

We carry out simulations of two discs with U oc and 

radial extent 1 < R < 100 AU. Since the aim of this pa¬ 
per is to investigate the radial movement of the disc mass 
following the initial fragmentation, we set up the discs such 
that they do f ragment. Disc s are more likely to fragm ent 
at large radii (|Rafikov gOOgj^ Clarkd [ 2 OO 9 I : iBolevI li^OPOl ) . if 


the disc mass is high (iRafikov 1200^1 . or if the metallicity is 
low (ICai et al ] l2006l : iMerii fc Bat,ell20fnl : iR.ogers Wadsle^ 
I 2 OI 2 L We assume the initial and boundary temperature pro¬ 
file to be T oc and set the absolute value of the tem¬ 

perature from the luminosity, L, of the central star (using 
L — where a is the Stefan-Boltzmann constant). 

The lumi nosity is determin ed using the stellar evolution 
models of ISiess et al.l (|2000l l for a low metallicity Z = 0.01 
at an age of 1 Myr (i.e. at a sufficiently early age when 
a disc would have been present around the central star). 
The resulting surface temperature of the discs modelled are 
~ 10 — 50 K (see Figure [1]), in a greement wit h previous sim¬ 
ilar fragmentation simulations: IBoIcv (l2009ll sets the disc’s 
i rradia tion temperature to 30K. IStamatellos &: WhitworthI 
(l2009h have one that var ies with radius and is 0(10) K, and 
NavakshinI (|201lh set their ambient gas temperature 
to lOK. It is important to note that while the luminosity 
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Triggered fragmentation in GI dises 3 


may be important for disc fragmentation (since it affects 
the disc energetics), this investigation is concerned with the 
dyn amics of the d i sc rath e r than its thermodynamics. 

iMeru Batel (|2011al . I 2 OI 2 I I have called into question 
the detailed thermodynamical conditions under which frag¬ 
mentation can occur. However iLodato Clark3 (|201lh sug¬ 
gest that the results concerning the dynamics of such discs 
are more likely to be robust since the resolution that is re¬ 
quired to accurately model such large-scale processes are not 
as stringent as that required to resolve fragmentation. Since 
we are more concerned with the dynamics of the discs, i.e. 
how the mass in the disc moves in response to the presence 
of a fragment and the potential for further fragmentation 
rather than the actual process of fragmentation itself, we 
simulate these discs using 250,000 particles. 

We firstly simulate a I.IM© disc surrounding a l.SM© 
star with luminosity L = d.SL©. We model this disc with 
opacity values that are 0.9 times the interstellar Rosseland 
mean values. We also present the results of a second simu¬ 
lation with a disc mass of O.OM© around a l.OM© star with 
L — 2.4L0. To ensure that the disc does initially fragment, 
we simulate it with a low opacity (0.5 times the interstellar 
Rosseland mean values). The setup is such that the min¬ 
imum value of the stability parameter (at the outer disc 
edge) is Qmin = 0.75 and 0.7 in Simulations 1 and 2, re¬ 
spectively. However, the discs do heat up due to work done 
and viscous heating and hence the stability parameter does 
increase from the initial values. 


4 RESULTS 

Figure!^ shows the radial velocity and surface mass density 
rendered images of the disc in Simulation 1 prior to the for¬ 
mation of any fragments when the gravitational instability 
has developed. The radial velocity plots are determined by 
projecting all the SPH particles onto a single x-y plane. The 
left panel shows that the movement of mass in the disc is 
both inwards and outwards but that the radial velocity of 
the gas is not extremely large in either direction. 

Figure [3] shows the radial velocity and surface mass den¬ 
sity rendered images of the disc a short while after the first 
fragment forms (the fragment forms at ~ 54.5 AU at a time 
t ~ 1.30 ORPs). This figure shows that mass is moving 
both outwards and inwards in the vicinity of the first frag¬ 
ment (the inwards motion is represented by the black region 
in the left panel). The magnitude of the radial velocities are 
up to a factor of ~ 10 larger than before fragmentation (com¬ 
pare solid and dotted lines in Figure [U and hence the disc 
is much more dynamic with more extreme mass movement. 
This movement into the inner regions causes mass to accu¬ 
mulate onto an inner spiral arm (indicated by an arrow in 
the right panel of Figure [3]). This spiral arm then fragments 
a short while later (at ~ 24 AU at a time t — 1.76 ORPs; 
Figure O. 

Figure [6] shows the azimuthally averaged stability pro¬ 
file for the disc just before the formation of the first and 
second fragments, as well as a 2D map of the stability pa¬ 
rameter in the disc before the first and second fragments 
form. The 2D stability values are calculated using the mass 
and internal energy of the SPH particles that are closest to 
the x-y position plotted, i.e. by projecting the SPH parti- 



Radius [AU] 

Figure 4. Azimuthally averaged radial velocity profile in Simu¬ 
lation 1 shortly before the formation of the first (t = 1.19 ORPs; 
solid line) and second (t = 1.75 ORPs; dotted line) fragments. 
The radial velocities increase by up to an order of magnitude 
after the first fragment forms compared to before. The dashed 
line shows the radial velocity of the disc at t = 1.75 ORPs but 
in the disc where the formation of the first fragment has been 
artificially suppressed fSection 14.11) . Without the presence of the 
first fragment the movement of material in the disc is no longer 
enhanced. 

cles onto a single x-y plane. This is sufficient as firstly, the 
temperature of most of the particles lies within a factor of 
a few of the midplane temperature, especially in the outer 
disc which is close to being vertically isothermal (see Fig¬ 
ure [J) . Secondly, the stability parameter only depends on 
the square root of the temperature. Prior to the first frag¬ 
ment forming, regions of the outer disc are gravitationally 
unstable (the stability parameter, Q < 1) but at R < 50 AU, 
Q > 1, and thus is not so susceptible to fragmentation. How¬ 
ever, following the formation of the first fragment, the inner 
spiral becomes sufficiently dense such that its stability pa¬ 
rameter decreases, allowing it to be pushed (from a state of 
marginal stability where Q > 1) into a fragmentation state 
(at 24 AU). 


4.1 Would the inner disc have fragmented 
anyway? 

In the above section we suggest that the formation of the sec¬ 
ond fragment was a consequence of the presence of the first 
fragment and the resulting movement of disc gas. However, 
to show that the second fragment would not have formed 
regardless of the presence of the first fragment, we restart 
the simulation at an earlier time, t — 1.04 ORPs (shortly be¬ 
fore the increase in surface mass density begins at the radius 
where the first fragment forms), but remove all the particles 
that are within a 1 an annulus of where the first fragment 
forms, i.e. the particles within 53.5 < R < 55.5 an (totalling 
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Figure 2. Radial velocity (left panel) and surface mass density (right panel) rendered images of Simulation 1 at a time {t = 1.07 ORPs), 
shortly before the formation of the first fragment (at t = 1.30 ORPs). At this time the gravitational instability has developed and 
the radial velocity is neither at one extreme nor the other. Negative and positive values of the velocity indicate inwards and outwards 
movement, respectively. The central star is given by blue and black dots in the left and right panels, respectively. 


Z) 

< 



Figure 3. Radial velocity (left panel) and surface mass density (right panel) rendered images of Simulation 1 at a time, t = 1.75 ORPs, 
shortly before the second fragment forms. After the formation of the first fragment material moves radially inwards (marked by an arrow 
in the left panel), causing the inner spiral (indicated by an arrow in the right panel) to become sufficiently dense that it can fragment 
(Figure [5j. The central star and fragment are given by blue and black dots in the left and right panels, respectively. 


O.O19M0 which is 1.8% of the SPH particles at that time). 
This artificially stops the first fragment forming and allows 
us to see if the second fragment would have formed anyway. 
We run the simulation until 8.9 ORPs (i.e. 5 times longer 
than it took until the second fragment formed in the orig¬ 
inal simulation) and find that the inner fragment does not 
form at all. In addition, the radial velocity rendered images 
show a quiescent disc, i.e. no mass movement takes place, 
similar to what is seen before any fragmentation occurs in 
the original simulation (Figure [2|). 

To ensure that this artificial removal of particles from 
the annulus does not cause the gas from the inner disc to re¬ 
plenish this artificial gap (and thus artificially decrease the 
surface mass density of gas near to where the inner frag¬ 
ment forms) we observe the stability parameter plot of the 
inner disc. Figure [7| shows the stability profile of this disc 
at various times (before the particles are removed and af¬ 


terwards) and shows that the stability profile at 24 an 
marginally decreases over time. If the artificial removal of 
particles caused material to move away from the disc region 
at 24 an the surface density would decrease and the sta¬ 
bility parameter would increase. The fact that this is not 
happening shows that the lack of fragmentation cannot be 
due to this. This provides further confidence that the reason 
for the second fragment forming is due to the presence of the 
first fragment. Furthermore, Figure [8] shows the radial veloc¬ 
ity rendered plot of the image at the end of the simulation 
(at t = 8.9 ORPs) and shows that the mass movement in 
the disc is as quiescent in the inner regions as it was before 
any fragmentation occurred (Figure [2l left panel). 
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Figure 5. Radial velocity (left panel) and surface mass density (right panel) rendered image of Simulation 1 at the time when the second 
fragment forms {t = 1.76 ORPs). This fragment forms out of a spiral arm that becomes dense due to the movement of mass radially 
inwards following the formation of the first fragment (Figure [S]). The central star and fragments are given by blue and black dots in the 
left and right panels, respectively. Following the formation of the second fragment the disc becomes much more dynamic. The presence 
of multiple fragments causes several regions to have radial velocities that are higher than they were prior to fragmentation (c.f. Figure [2] 
left panel and Figure [4|). 


4.2 Simulation 2 

Simulation 2 shows a similar behaviour to that of Simu¬ 
lation 1. The disc becomes unstable but is quiescent before 
any fragmentation occurs (Figure [9l top row). The disc frag¬ 
ments first at ~ 33 an (at t = 1.94 ORPs; Figure [9] second 
row) after which it becomes more dynamic with more mass 
moving in both directions (Figure[9l third row). The inwards 
movement causes gas to pile onto an inner spiral (indicated 
by an arrow in Figure and fragment (at R ~ 10 an at 
t — 2.2 ORP; Figure [9] bottom row). Figure fTOl shows that 
before the first fragment forms the inner disc is stable and 
is not expected to fragment (Q > 1) but at a late time this 
region is pushed into a state of instability. 

As with Simulation 1, to test whether the inner frag¬ 
ment would have formed anyway we restart the simulation 
at t = 1.3 ORPs but remove the particles from the annulus 
between 32 < R < 34 an (totalling O.O 2 IM 0 which is 2.6% 
of the SPH particles at that time) so that the formation of 
the outer fragment is suppressed. We run the simulations 
until 11.1 orbits (more than 5 times the time it took until 
the second fragment formed in the original simulation) and 
find that no fragments form. Figure ITT] (left panel) shows the 
stability map at the end of the simulation and shows that 
the disc is stable, and the radial velocity rendered image of 
the disc shows it to be calm and quiescent (Figure [TTl right 
panel), as it was before any fragmentation occurred in the 
original simulation (Figure [9l top left). 


5 DISCUSSION 

Disc fragmentation is known to occur at large radii in young 
discs. However the subsequent evolution of a disc immedi¬ 
ately after fragmentation and the potential for further frag¬ 
mentation as a result of a fragment forming is typically un¬ 
derstudied when considering this mechanism. In this paper 
we show that the movement of mass as a consequence of 


disc fragmentation at large radii can potentially cause an in¬ 
ner spiral in a disc to become sufficiently dense that it may 
fragment. The interaction between the disc gas and the frag¬ 
ment can be thought of as analogous to the theory behind 
gap-opening by planets. A region of gas just interior to the 
fragment transfers angular momentum to the fragment. Its 
magn itude is approximately given by (see e.g. iLubow Idal 

ionl): 


|AJ| 



( 2 ) 


where M is the mass of the region of gas being considered, 
Qf and Mf are the angular frequency and fragment mass, 
respectively, M* is the mass of the central star, a is the 
fragment’s radial location in the disc and x is the distance 
between the fragment and the region of gas being considered. 
This causes the region of gas to move inwards. Similarly 
angular momentum is also transferred from the fragment 
to the gas just exterior to it causing the gas to move to 
larger radii. The inwards motion potentially causes an inner 
spiral to become sufficiently dense such that it can become 
unstable and fragment. This is particularly likely if the inner 
spiral is at a marginally stable state (with Q 1) and just 
requires an additional push to cause it to become unstable 
and fragment. The inwards motion is analagous to a high 
mass accretion rate through the disc (iHavfield et al. 2011) 
or a high infall rate (IVorobyov Basull2005l : iKratter et al.l 
bnnsi: lBolevll200fll : IZhii et al.ll2012l i which increases the local 
disc surface density and decreases the local stability, causing 
frag mentation to oc c ur. 

I Vorobyov et al.l (|2013l i suggest that mass infall may 
cause periodic bursts of fragmentation. It is also possible 
that the infall causes the first fragment to form, following 
which further fragmentation is triggered (or accelerated) ac¬ 
cording to the mechanism in this paper. 

It is important to note that once multiple fragments 
form, the presence of all the fragments, as well as the sec- 
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Figure 6. Top panel: Azimuthally averaged stability profiles of the disc in Simulation 1 shortly before the first (solid line; t = 1.19 
ORPs) fragment forms and at two instances before the second fragment forms (at t = 1.62 and 1.76 ORPs; dotted and short dashed 
lines, respectively). Before any fragments form, part of the outer disc (beyond 50AU) is unstable and can fragment while the inner disc 
at 24 an (where the second fragment later forms) is stable {Q > 1). The mass movement into the inner parts of the disc following the 
formation of the first fragment causes the inner spiral to be pushed into a state of instability such that it fragments (at 24 AU). Note 
that the stability profile does not show that the disc is unstable at this location. However such azimuthally averaged plots wash out the 
local properties. Bottom panels: Local stability parameter values before the first and second fragments form (at t 1.19 and 1.76 ORPs, 
respectively). The black circles indicate the radius at 24 an. Before the first fragment forms the value of the stability parameter is > 1, 
and hence stable, inside 50AU (bottom left). Before the second fragment forms regions of the annulus at 24 an are definitely unstable 
with Q < 1 (bottom right). The second fragmentation is therefore unlikely to have occurred without the aid of the mass movement. The 
central star and fragment are given by black dots. 


ular evolution of the disc, will play a part in determining 
how the disc mass moves. Indeed, when run for longer, more 
than two fragments do form. Previous simulations suggest 
that a period of chaotic behav iour may occur when m ul- 
tiple objects co-exist in a disc (jPierens Nelsonll2008l ). It 
is beyond the scope of this paper to characterise the de¬ 
tailed mass movement that may occur in any one particular 
case since each disc is likely to be dynamically different to 
one with different disc properties or a slightly varying evo¬ 
lution of fragments. We stress that the focus here is that 
the mass movement in a disc is much more dynamic after 
fragmentation than before, and that this dynamical activity 
may cause subsequent fragmentation. We also stress that 


since the details into the fragmentation of s elf-gravitating 
discs h a ve recently been called into question (iMeru Bate 


2011a 

bl. 2 OI 2 I: Lodato & Clarkd I 2 OIII: IPaardekoooer et al. 

201 ll 

Rice et al.ll2014l) with further comolications when con- 


sidering that fragmentation is also somewhat a stochasti c 
process (|Paardekooper|[201^ : iHopkins Christiansen|[2013l ) , 
the emphasis here is not on fragmentation itself, but the 
dynamical movement of the disc mass following the initial 
fragmentation stage, which can then cause an inner spiral to 
become sufficiently dense such that it can potentially frag¬ 
ment. 

iMeru Batel (|2012l ) showed using their simulations of 
fragmenting self-gravitating discs using a simplified cooling 
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Figure 9. Radial velocity (left panels) and surface mass density (right panels) rendered images of Simulation 2 before any fragments 
form (at t = 0.92 ORPs; top row), when the first fragment forms (at t = 1.94 ORPs; second row), shortly before the second fragment 
forms (at t = 2.18 ORPs; third row) and after the second fragment forms (at t = 2.20 ORPs; bottom row). The disc is initially relatively 
quiescent. After the first fragment forms the disc becomes dynamic, causing gas to move more rapidly in both directions. The inwards 
movement results in mass piling onto an inner spiral causing it to fragment. The central star and fragment are given by blue and black 
dots in the left and right columns, respectively. 

© 0000 RAS, MNRAS 000 , 000-000 









































8 Farzana Meru 



Figure 10. Local stability parameter values of the inner disc in Simulation 2 just before the first (left panel) and second (right panel) 
fragments form (at t 1.9 and 2.2 ORPs, respectively). The black circle indicates the radius at 10 au. Before the first fragment forms the 
inner disc is stable and cannot fragment. However just before the second fragment forms regions of the annulus at 10 au are definitely 
unstable with Q < 1, allowing the disc to fragment. The second fragmentation is unlikely to have occurred without the aid of the mass 
movement since the stability parameter value is > 1 at 10 au at the earlier time. The central star and fragment are given by black 
dots. 
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Figure 11. Local stability parameter rendered (left panel) and radial velocity rendered (right panel) plot for the disc at the end of 
Simulation 2 (at t = 11.1 ORPs) where an annulus of particles has been removed. The inner disc is stable (Q > 1) and does not fragment if 
the outer fragment doesn’t form. In addition the disc is quiescent without significant mass movement, as it was in the original simulation 
before any fragments formed (Figure O top left) 


method that the values of the artificial viscosity parame¬ 
ters used in simulations can result in counterintuitively high 
dissipation, which may cause fragmentation to be inhibited. 
While it remains to be tested exactly how the dissipation 
affects the fragmentation in discs modelled with radiative 
transfer, we stress that the inhibition of fragmentation only 
means that our results are erring on the side of caution 
rather than over-predicting fragmentation. 

In the simulations presented here, since we have mod¬ 
elled a disc more realistically using radiative transfer and by 
considering the effects of stellar irradiation on the surface 
temperature of the disc, the initial and subsequent fragmen¬ 
tations take place in the outer and inner regions, respec¬ 
tively. We point out that if a simulation is carrie d out us¬ 
ing a simplified parameterisation for the cooling (|Gammiel 
l 2 nnil: as done so in many previous simulations), the first 


fragments would form at small radii (since the absol ute cool¬ 
ing timescale is shortest there; iMeru Batdlioilbl ) and the 
subsequent fragmentation would then occur in the outer 
regions. While in this case, the subsequent fragmentation 
would be expected to occur anyway, the fragmentation may 
occur faster due to it being aided by the movement of mass 
in the disc (though this time in an outwards direction). Fur¬ 
thermore, in discs modelled using radiative transfer which 
may subsequently fragment anyway either interior or exte¬ 
rior to the first fragment irrespective of the gas dynamics, 
the mass movement may cause the disc to fragment faster 
by aiding the destabilisation process. 
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Figure 7. Azimuthally averaged stability profiles of the disc in 
Simulation 1 that is modified to remove an annulus of particles 
at various times: just before any particles are removed it = 1.04 
ORPs; solid line), at the time when the first fragment formed in 
the original simulation {t = 1.30 ORPs; dotted line), at the time 
when the second fragment formed in the original simulation (t = 
1.76 ORPs; short dashed line) and much later in the simulation 
{t = 3.0 ORPs; long dashed line). The stability parameter does 
not change significantly in the inner disc at 24 an. This suggests 
that the fragmentation in the outer disc may be the cause of the 
inner disc fragmenting in the original simulation. 
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Figure 8. Radial velocity rendered image of the disc in Simu¬ 
lation 1 that is modified to remove an annulus of particles. The 
image is at the end of the simulation when no fragments have 
formed. The disc is quiescent in the inner regions, similar to the 
original simulation before fragmentation occurs (Figure [21 left 
panel). The disc is also much more quiescent in the inner disc 
compared to the original simulation when fragmentation occurs 
(Figure (3) left panel). 



5.1 In the context of planet formation 

Our results may have important implications for the theories 
of giant planet formation: whi le a disc may only initially 
fragment in the outer regions (jRafikovI |200^ : IClar ^ 120091: 
iBolevI EoOflh ■ planets may well form in the inner parts via 
gravitational instability provided fragmentation has already 
occurred in the outer disc first. 

Core accretion has historically been thought to oc¬ 
cur within a disc’s lifetime at small radii (up to 
5 — 10 AU). Gravitational instability, on the other 
hand, is tho u ght to occu r at large radii (> 50 AU; 

iRafikovI l2009l : Iciarkd l2009l : iMeru Batel 1201^ . Though 


a number of factors may influence these radii, e.g. disc 
mass, stel lar mass, meta ll icity (or equivalently, disc opac- 


ity) (e.g. ISafronov 19691: Whitworth & StamatellosI 

2006; 

Cai et al.l I 2 OO 6 I: ICossins et ahl 2010l: iMeru & Bate! 

2010; 

Rogers & Wadslevll2012l). an intermediate region aoDcars to 


exist (~ 10 — 50 AU) in which in situ formation by either 
mechanism is somewhat harder. R ecent advances in growth 
models at various size scales (jha mbrechts fc Jo hansenI 

I 2 OI 2 I : ICaraud et al.ll201^ : iHelled Bodenheimerll2014l i have 

helped to improve the possibility of in situ growth by core 
accretion at larger distances and the results presented in 
this paper show that subsequent planet formation by gravi¬ 
tational instability may well be a mechanism that may op¬ 
erate to form planets at smaller radii than expected. Thus 
both formation models have made progress in the interme¬ 
diate radial range where no one single formation method has 
historically been thought to dominate. 

lArmitage HansenI (| 19991 ) first presented the idea that 
the presence of a planet in a self-gravitating disc may trig¬ 
ger the formation of a fragment in another part of the disc 
(both interior and exterior to it) where the stability param- 
^ eter is 1 — 2. However, these discs were modelled using an 
■g isothermal equation of state (which favours fragmentation) 
— and by inserting a planet artificially. We show that triggered 
■■g fragmentation occurs even when considering more realistic 
■5 energetics and when the formation of the first fragment is 
modelled self-consist ently. Our results are never theless Con¬ 
'S sistent with those of lArmitage &: HansenI (| 19991 ). 

We also note that the simulations presented here are 
of high mass discs. The mass of the fragments when they 
first form range between 0.7 — 1.2Mj. However, by the time 
the second fragments form the first fragments have grown 
to ~ 62Mj in both simulations. While we expect that mass 
movement would also take place in lower mass discs, the ex¬ 
tent to which the inner spirals become sufficiently dense will 
be different. Equation [2] shows that the amount of angular 
momentum transfer depends not only on the fragment mass 
but also on the mass of the region of gas being considered 
(both will ultimately depend on the disc mass). Therefore, 
we expect that the mass movement in a lower mass disc will 
be less violent. While we have shown the mass movement 
to be important in high mass discs which are in their very 
young phase, we present this as a possible mechanism by 
which further fragmentation at smaller radii may occur after 
an initial fragmentation stage in lower mass self-gravitating 
discs, where the disc mass is of the order of 10% of the stel¬ 
lar mass. However, this needs to be verified with numerical 
simulations. Furthermore, in the context of planet formation 
it is possible that the fragments could lose mass to become 
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planetary mass objects (iBolev et al.ll2010l : lNavakshinll201Clh . 
However, our simulations use sink particles which can only 
grow with time so cannot model this in detail. 


6 CONCLUSIONS 

We carry out three-dimensional radiation hydrodynamical 
simulations of self-gravitating discs to investigate the dy¬ 
namics of the disc following the initial formation of a frag¬ 
ment. In particular, we self-consistently explore the forma¬ 
tion of the first fragment and the resulting potential for sub¬ 
sequent fragmentation. We find that the movement of disc 
mass after the first fragment forms is much more dynamic 
than prior to fragmentation as the gas radial velocity in 
parts of the disc can increase by up to a factor of ~ 10 in 
both the inward and outward directions. 

The inwards movement of mass can cause an inner spi¬ 
ral arm in a gravitationally unstable disc to become suffi¬ 
ciently dense, pushing it from a stable or marginally stable 
state into a state of instability such that it may fragment. 
This may potentially form fragments at radii smaller than 
typically expected in protoplanet ary discs (^50 AU). This 
subsequent fragmentation may provide a possible method 
by which planets may form in the intermediate radial region 
10 — 50 AU) where no one in situ formation method has 
historically dominated. 
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